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, We present and discuss the derivation of a nonlinear non-local integro-difFerential 

Ci ' equation for the macroscopic time evolution of the conserved order parameter p{r, t) 

W), of a binary alloy undergoing phase segregation. Our model is a d-dimensional lattice 

^ ^ gas evolving via Kawasaki exchange dynamics, i.e. a (Poisson) nearest-neighbor ex- 

, change process, reversible with respect to the Gibbs measure for a Hamiltonian which 

^ ' includes both short range (local) and long range (nonlocal) interactions. The nonlo- 

^ , cal part is given by a pair potential 7"^ J(7|a; — y|), 7 > 0, a; and y in Z*^, in the limit 

7 — > . The macroscopic evolution is observed on the spatial scale 7^^ and time 
scale 7^^^, i.e., the density, p(r, i), is the empirical average of the occupation numbers 
' over a small macroscopic volume element centered at r = 7a;. A rigorous derivation 

is presented in the case in which there is no local interaction. In a subsequent paper 
(part II), we discuss the phase segregation phenomena in the model. In particular 
we argue that the phase boundary evolutions, arising as sharp interface limits of the 
family of equations derived in this paper, are the same as the ones obtained from the 
corresponding limits for the Cahn-Hilliard equation. 
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1. Introduction 

The process of phase segregation, following a quench (sudden cooling) of a sys- 
tem from high temperature, where the system has a unique uniform equilibrium 
phase, into the miscibility gap where two (or more) phases can coexist, is variously 
known as spinodal decomposition, nucleation, coarsening, etc. . . It concerns the 
tendency of the different phases to segregate, creating larger and larger domains 
of approximately homogeneous single phase regions. The problem is of great prac- 
tical importance in the manufacturing of alloys where the degree of segregation 



influences the properties of the material. The mathematical description of the time 
evolution of the local macroscopic order parameter in such systems, e.g. the differ- 
ence in the concentration of A- and B-atoms in a binary A-B alloy, is commonly 
given by nonlinear fourth order equations of the Cahn-Hilliard type [4]. 

These equations appear to capture much of the phenomena. In particular, their 
numerical solutions show good agreement with experiments and with computer 
simulations of the Ising model (thought of as a binary alloy) evolving via Kawasaki 
exchange dynamics. This agreement relates both to the appearance and shape of 
segregated domains, which seem to exhibit a self-similar structure, and also to the 
quantitative behavior of the characteristic length describing these structures, which 
seems to grow like t^^^, where t is the time. While this agreement is certainly satisfy- 
ing, the original Cahn-Hilliard equation (CHE), or the modifications of it proposed 
so far [9], do not seem to to arise as exact macroscopic description of microscopic 
models of interacting particles, such as the Ising model with Kawasaki dynamics 
(the CHE can however be derived from certain mesoscopic Ginzburg-Landau con- 
tinuous spin models [2]). This is unlike some other physically motivated equations, 
e.g. the diffusion equation and the Boltzmann equation, which can be derived from 
idealized microscopic models in suitable limits [7,29]. Such derivations are both of 
intrinsic interest and also indicate something about the range of applicability of 
the macroscopic equations. The latter might be particularly relevant for the CHE, 
where all that is known mathematically about the behavior of the solutions is re- 
stricted to the late stages of the coarsening process when the evolution is assumed 
to be dominated by the motion of sharp interfaces between well formed domains of 
the pure phases. It is far from clear how much this singles out the CHE from other 
possible equations describing phase segregation. 

In this note we rigorously derive a macroscopic equation, describing phase segre- 
gation in microscopic model systems with long range interactions evolving according 
to stochastic Kawasaki dynamics with nearest neighbor exchanges. We will then, 
in Part II [14], study the interface motion obtained from the derived macroscopic 
equation, in the sliarp interface limit, by means of formal matched asymptotic ex- 
pansions of the solution of the macroscopic evolution equation (see e.g. [3,5,18,24]). 
By sliarp interface we mean the limit in which the phase domains are very large 
compared to the size of the interfacial region, i.e. denoting by L the typical size 
of the domains, we look for results in the limit L — > cxd. The time will have to be 
properly scaled as well, typically as some integer power of L, according to the type 
of initial condition and the choice of the temperature. Our conclusion there is that, 
from the sharp interface viewpoint, the equation derived from a particle system and 
the Cahn-Hilliard equation are essentially equivalent. 

The models we consider are dynamical versions of lattice gases interacting via 
long range Kac potentials, also known as local mean field interactions. The equilib- 
rium properties of these systems are well known ([16,21,25]). They provide micro- 
scopic models in which the van der Waals or mean field description of phase transi- 
tion phenomena, which is in good qualitative (or even quantitative) agreement with 
experiments away from the critical point, holds exactly. This includes metastabil- 
ity phenomena [25] . The corresponding dynamical models which we study here are 
sometimes called local mean field Kawasaki dynamics. They can be described in 
words as follows: each particle hops (at a random time) from a site of the lattice 
Ia'^ to one of its 2d neighboring sites with a rate which depends on the particle con- 
figuration. These rates are chosen to satisfy detailed balance (reversibility) with 
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respect to the Gibbs measure having the specified interaction between the particles 
[29]. 

In the simplest model we consider here, there is only a long range (Kac type) 
potential. We will also discuss, but not investigate in detail, the case in which ad- 
ditional short range interactions are also present. Further work on the same model 
and related ones can be found in [12,20,32]: these papers focus on the diffusive 
regime, i.e. on the region in which there is only one phase, but versions of the 
integral equation on which we are focusing are already present there. The case 
without a conservation law, Glauber dynamics, has received much attention, see [6] 
and references therein. 

The precise definition of the model is given in the next section. We show there 
that the evolution of the macroscopic density (in the simplest model) is given in 
terms of a second order integro-differential equation (2.16). We then show that this 
equation can be written in terms of the gradient flux associated with the classical 
local mean field free energy functional ([21,25]) and a density dependent mobility. 
This allows us to make a direct connection between the properties of the solution of 
the derived evolution equation and the equilibrium phase diagram as well as with 
the solution of the CHE. In section 3 we argue that the gradient structure for the 
macroscopic evolution law should hold generally in systems with long range Kac 
potentials and arbitrary additional short range interactions. In Section 4 we make 
a remark on the case in which a weak external field is present. The proof for the 
simple case is given in section 5. 

2. The particle model and its hydro dynamic limit 

The particles live on the (i-dimcnsional lattice A-y = {1, 2, . . . , [7"''^]}'', where 
7 > is a small parameter and [r] denotes the integer part of the real number r. 
We impose periodic boundary conditions on A^. Each site of A^ is cither occupied 
(1) or empty (0), hence a particle configuration is an element ry of fi-y = {0, 1}^^ 
and the latter is endowed with the product topology. The dynamics is specified by 
giving an initial condition rjo G fi-y, which may be random, i.e. a measure on the 
Borel sets of 0-y, and some (stochastic) evolution rules which will define r]t for any 
t positive. Our aim is to have a Markovian dynamics for which the Gibbs measure 
associated with a given Hamiltonian and a given total particle number is the unique 
reversible time invariant measure at a fixed temperature. The Hamiltonian is a real 
valued function defined on and we take it to be the sum of two terms: 

H = Hs + H^ (2.1) 
Hs = -l Yl K{x-y)r^{x)r^{y) (2.2) 

^''Ai{^-y)HxUy) (2.3) 

in which J is a smooth function (C°°) from the d-dimensional unit torus T'^ to 
the real numbers, such that J(r) = J{—r), and K{x) = if > R, for some R 
independent of 7. The term in (2.3) will be called nonlocal, while the one in (2.2) 
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will be called local or short range. The Gibbs measure with Hamiltonian at the 
temperature (/3 > 0) and total number of particle A?" e Z+ is defined as 



where 

Z^(iV)= Yl exp(-/3i/(r7)). (2.5) 

(2.4) is a probability measure over O:^ = {77 G : XlxeA ^/(^) = The 
stochastic process {'r]t}t>o is the Poisson jump process [22,29] generated by the 
operator 

Lyf{v)= E (2-6) 

where / is a real valued (bounded) function on Qj, 

{r]{x) \iz = y 

r}{y) iiz = x (2.7) 
r]{z) otherwise 

and 

c^{x, y; r/) = $ {(3 [Hirj'^'y) - Hirj)]} (2.8) 

if |a; — y| = 1 and is zero otherwise. Here $ : R — > R"'" is twice differentiable in a 
neighborhood of and satisfies the detailed balance or reversibility condition, 

^E) = exp{-E)^{-E) (2.9) 

for all E E M.. We also assume $(0) = 1: the case $(0) G (0, 00) can be recovered 
with a time change. 

Loosely speaking, the process t^t can be described in the following way: if at time 
t the configuration is r]t , the probability that in the time interval [t, t + At] the sites 
x,y {\x — y\ = 1) exchange their occupation numbers is 

c^{x,y;r]t)At + 0{{At)^). (2.10) 

We note that if r]{x) = rj{y), then an exchange between x and y does not modify 
the configuration rj and it is thus possible to interpret the dynamics in terms of 
particles which attempt to jump from x to but the jump is performed only if 
the site y is empty: this type of dynamics is said to have an exclusion rule, that 
is the particles have an on site hard-core repulsion. A detailed construction of 
this process in terms of Poisson jump processes is given on page 158 of [29]. The 
configuration space we are working on (Q^) is finite and this avoids the difficulties 
connected to defining such a dynamics on an infinite state space (Chapter 1 of [22]). 
In Section 1.2 (part II) of [29] it is shown that for / and g, bounded functions on 

j 9iv)L^fiv) d/x?(r/) = j fiv)L^9iv) d//?(r/) (2.11) 



This property is called reversibility and it is a direct consequence of (2.9). In 
particular (2.11) implies that 



J L^firj) d//?(r7) = 



(2.12) 



We will use the following notation: the generic initial condition is a probability- 
measure ji on Q^. The law of the process {r]t}t>o with initial condition /i will be 
denoted by (-Ej^ for the expectation) . The process PJ^ is linked to the semigroup 
generated by via the formula Ei^{f{r]t)) = f{exp{L^t)f){ri)iJ,{dr]). Equation 
(2.12) implies that for any t >0 and any / bounded 



dt ^ 



(2.13) 



which means that fj,^ is invariant under the process generated by L^. 

The hydrodynamic limit. We are interested in initial states /j,^ such that, when 
7 — > 0, //^ resembles more and more a profile po, where po is a measurable function 
from the (i-dimensional unit terms T*^ to [0, 1], stretched by 7~^. More precisely, 
we say that {P'-yj-yyo is an initial condition associated with po if, for any continuous 
function from to M and every 5 > 0, 



lim n~ 

7^0 ' 



7 



y (l){'yx)r]{x) — / (j){r)po{r)d7 



> 5 



(2.14) 



The condition (2.14) is clearly satisfied if p,^ is such that {r]{x)) = r]{x)dp^{r]) 
= poi^x) for all X in and the occupation number of the sites are independent. 
On the other hand, the initial condition concentrated on the chessboard configu- 
ration {r]{x) = 1 if the sum of the components of x is even, r]{x) = otherwise) 
is also obviously associated with po = 1/2 and many other examples can be easily 
constructed. Loosely speaking, when we say that the particle system has a hydro- 
dynamic limit, we mean that (2.14) holds also at later times, if we replace po by the 
solution of a suitable hydrodynamic evolution equation (in our case we will have an 
integro-differential equation), with initial condition po. In Section 4 we prove the 
following Theorem: 

Theorem 1. The Hydrodynamic limit without short range interactions. Set K{x) 
= so that the Hamiltonian coincides with H^. Let p^ he an initial condition 



associated with po G . 
function (p from T'^ to M 



Then, for any t positive , S > 0, and any continuous 



lim 

7^0 ~ 



7 



V (l){-fx)r]tj-2{x) - / (j){r)p{r,t)dr 



>5 \ =Q 



(2.15) 



when p{r, t), r E T^ and t e [0, oo), is the unique classical solution of the equation 
dtp{r, t) = V • [Vp(r, t) - (3p{r, t){l - p{r, t)) J^, VJ(r - r')p{r')dr'] 



(2.16) 



p(r, 0) = po(r) 



It is now an observation, which at first sight appears surprising, that (2.16) can 
be rewritten in the form 



dp{r,t) 
dt 



5p 



(2.17) 



where 



^\P) = --J s{p{T))dT -Iff J{t- r')p{T)p{r')dTdT' (2.18) 



with 



and 



s{p) = -plogp - (1 - p) log(l - p). 



a\p) = f3p{l-p) 



(2.19) 



(2.20) 



Rewriting (2.16) in the form (2.17), is recognized as the free energy functional 
and cr^{p) as the mobihty of our model without the long range interactions. This 
allows us to connect our equation with that of the CHE. Before doing that we 
rewrite T^{p)^ up to an irrelevant additive constant, in the form 

J^^{p) = f fl!{p{r))dr + 1 f J{r- r') [p{r) - p{r')f drdr' (2.21) 



where 



/c(P) 



J(0) 



P- 



+ (plogp +(l-p)log(l-p)) (2.22) 



and J(0) = / J(r)dr. 

Note that if p is constant the second term in (2.21) vanishes, so that /° is the 

constrained equilibrium free energy density of a homogeneous system see [21, 25]). 
fcip) is in fact the correct equilibrium free energy density as long as /3 < /?c = 
1/Tc = 4/J(0). If /3 > /3c, fcip) has a double minimum at p — p^, the two 

nontrivial solutions of log(p/ (1 — p)) = (3J{0){p — 1/2), and the correct free energy 
is then obtained by the double tangent construction. A heuristic derivation of (2.16) 
and an explanation of (2.17) is given in Section 3, where the generalization to the 
case K ^ is considered. (2.17) has the same structure as the CHE, whose various 
forms correpond to different c^ip) and (p) chosen by different authors [4,5,18,26]. 
What is common to the different CHE's in the literature is that the second term 
on the right side of (2.21) is of the form |C /(Vp)^dr where C > is related to 
the surface tension. This can be thought of as expanding the term in (2.21) and 
keeping only some terms, which is reasonable when the scale on which p varies is 
large compared to 7"^. We shall discuss the relationship between the solutions of 
(2.16) and the CHE in [14], see also [13]. 



3. The general case: K ^ 

We start by giving a heuristic explanation of the result in Theorem 1, which 
is in fact a sketch of its proof. We are for the moment still in the case K = 0: 
if P = the particle system reduces to the symmetric simple exclusion process 
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(SEP), i.e. all the particles are performing exchanges with rate one, so their only- 
interaction is given by the exclusion rule. As it is straightforward to verify, the 
Bernoulli measures fip with uniform density p G [0,1], under which the random 
variables {'n{x)}xeA^ are independent and J ?7(0)d/Xp(?7) = p, are invariant for the 
SEP dynamics and the hydrodynamic limit for the SEP is simply given by the heat 
equation (see e.g. [17]), as in (2.16) with /? = 0. Moreover, we observe that in the 
case /? > 0, 



H(rj-'-+<^) - H{rj) = 7 {r}{x + e) - rj{x)) 



7''^r/(^)(e-VJ)(7(a;-^)) 



+ 0(7^) 



= 0(7) (3.1) 

for all a; e and e a unit vector in ((3.1) is derived in Lemma 2, Section 5), 
so that 

c^(x, X + e; 77) = 1 - I [H{r^-^-^^) - H{r^)\ + 0(7^) 

= 1 + 0(7) (3.2) 

where we used (3.1) and the fact that (2.9) and $(0) = 1 imply $'(0) = -1/2 . 

Formula (3.2) indicates clearly that the dynamics with /3 > is a weak perturba- 
tion of the /? = dynamics. In particular the dominant SEP dynamics will enforce 
at time ^7"^ [t > 0) local equilibrium with respect to its invariant measure, i.e. the 
state of the system will locally (on spatial scales shorter that 7~^) be very close to 
the Bernoulli measure pp and p will be varying on the macroscopic scale 7"^. The 
perturbation term in (3.1) generates a force term given by the negative gradient of 
the energy density at r, 

F(r) = V J(r - r')p{r')dr'^ (3.3) 

This gives an extra contribution to the macroscopic current equal to this force times 
the mobility cr° = /?p(l — p). In this expression /? measures the intensity of the bias 
in the exchange rates and p{l — p) gives the rate at which the exchanges actually 
take place when a particular site is chosen at random, since the system is locally 
described by the Bernoulli measure pp. This explains the form of (2.16). 

That (2.16) can be transformed into (2.17) is clearly due to the fact that s{p) in 
(2.19) satisfies the relation, 

-^s"ip)a\p) = 1 (3.4) 

for all p G (0, 1). This is no accident but, as will now show, (3.4) follows from the 
Einstein relation between fluxes and forces for this system. 

This will become clearer if we consider the general case K ^ 0. We call the 
system with only the local Hamiltonian Hg the reference system. The equilibrium 
free energy f^qip) associated with this reference system at an average density p G 
[0, 1], which depends also on /3, is uniquely defined in the thermodynamic limit (see 
[27, 28]). We shall further assume that feq{p) is strictly convex and real analytic, 
which implies that there is no phase transition for the equilibrium reference system, 
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associated with Hg at the temperature 1/p. Morever we are assuming that for each 
p G [0, 1] there exists a unique, translation invariant and ergodic, infinite volume 
limit Gibbs state fXp such that / r]{0)dfXp{rj) = p. All these properties, which are 
to a certain extent equivalent, are known to hold if (3 is sufficiently small (see [27, 
28]). 

We now claim that the correct macroscopic evolution law for the case K ^ 
should be 

dt 

where (Ts{p) is the mobility of the reference system. 

Hp) = J /c(p(r))dr J J{r- r) [p(r) - p(r')]' drd/ (3.6) 

and we have defined the constrained equilibrium free energy density (in analogy 
with [21]) as 

/c(p) = -^(p-0' + /|,(p) (3.7) 

which is clearly a straightforward generalization of (2.17). 

The definition of ag for a system with short range interactions in the one phase 
system can be found in Part II of [29], formulas (2.27), (2.71) and (2.72) (in our 
notation formulas (2.27) and (2.71) must be multiplied by /?). It is given by a Green- 
Kubo formula and is related to the diffusion coefiicient of the reference system Ds (p) 
by the relation, 

<Js = XsDs (3.8) 

where Xs{p) us the inverse of the derivative of the chemical potential, i.e. the com- 
pressibility of our reference system is defined as 

X. ^ T^ir^ = ^ E ( / V{XU0W;{V) - P') (3.9) 

While a complete proof of (3.7) for general short range interactions is lacking at 
the moment, some particular cases can be handled [1]. In order to justify (3.5) on 
a heuristic level we will draw an analogy with Section 3.4 of Part II of [29]. There 
a linear response argument is developed for the system with Hamiltonian 

Hviv) = + Yl y{^x)v{x) (3.10) 

where V{r) is a smooth function from T*^ to M and e > is a small parameter. This 
is the Hamiltonian of a weakly driven lattice gas: the dynamics is defined by the 
Poisson rates Ce,v{x,y',r]) — ^{(3Ax,yHv{r])) for \x — y\ — 1 (and Ce,v{x,y;r]) — 
otherwise). Hence we have 

c^^vix,x + e;r]) = (3.11) 

$(/3A,,,+ei?.(?7)) - /3e$'(/3A,,,+eifs(r7))(r7(x + e) - v{x)){e ■ VV){ex) + 0{e^) 

In the hydrodynamical scaling limit, with r scaled as ex, e — > 0, external force will 
vanish like e. Hence it can be argued, see in the macroscopic continuity equation 
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(3.5) 



for the density dp/dt = —V • J, the current J should be the sum of a diffusive 
term, unaffected by the weak external force, and a drift term proportional to F = 
-VV{r), 



J=-Ds{p)Vp + a{p)F (3.12) 

with cr(p) a mobility to be determined. To obtain a we observe that the Gibbs 
measure /ly associated with Hv{r]) at the temperature l/j3 is a stationary state 
of our system. This implies that dp{r,t)/dt = whenever p(r, 0) = Peq{f) is the 
density obtained from jjiy. In the limit e — > 0, peq is the density which minimizes 
the functional [21] 

^v{p) = I [A(p(r)) + V{r)p{r)] dr (3.13) 

with a total mass constraint. In particular, at the minimum, dJ-y/Sp = const. But, 
according to (3.12), p{r,t) satisfies the equation 

^=V-[L»,(p)Vp + a(p)VF] (3.14) 

so the requirement that dp/dt = at p{r) — Peq{f) then requires that a — DgX- 
This gives the described Einstein relation, showing that a = ag. 

In the model we are considering e = 7 and the external potential V is replaced 
by an internally generated one, — 7^^ XlyeA •^ili^ ~ y))viy)- Taking now the hy- 
drodynamic scaling limit, with the scaling parameter e set equal to 7 and letting 
7 — > 0, the previous analyses should remain valid with V{r) — J J{r — r')p{r')dr' 
in (3.14) which then gives (3.5). 

We note that if we look at a hydrodynamical space scale which is much larger 
than 7~^, say e = 7^, 5 > 1, then we get a purely diffusive equation for p{r, t). This 
is obtained from (3.5) by dropping the second term in (3.6). This equation holds 
rigorously outside the spinodal region, where the new diffusion constant becomes 
negative, see [20] and [12]. 

4. A REMARK ON THE EFFECT OF AN EXTERNAL DRIVING FORCE 

As one may expect, domain coarsening phenomena become richer and more 
complicated in the presence of an external driving force and some important points 
are still controversial (see e.g. [19,33]). One of the interesting problems is under- 
standing the scaling behavior of the cluster size; different behaviors are expected 
depending on whether the characteristic length is measured in the direction per- 
pendicular or parallel to the field. 

Here we simply observe that it is straightforward to extend the results of this 
paper to the case of a system on which a weak external force is acting. This system 
has a Hamiltonian 

H{ii) + yMvix) (4.1) 

where V is a smooth function from to R, as in Section 3. Remember that 
H = Hg + Hj (formula (2.2)). The discussion in Section 3 easily gives the following 
guess for the hydrodynamic limit of this system 

dp . 



dt 



f5J^ 
dV — + V 
\5p 



(4.2) 
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alternatively one can redefine J-{p) as J-{p) — J V{r)p{r)dr and the evolution equa- 
tion would still be (3.5). We observe that the proofs in Section 4 directly extend 
to this case, if = 0. 

The situation is similar in the case of an external force which is not the gradient 
of a potential F, e.g. a constant force E e W''. In this case the rates of the process 
would be 

$ iP{H{rj-'y) - H{rj) -^E-{x- y){rj{x) - rj{y))) (4.3) 
if \x — y\ = 1 and zero otherwise and the macroscopic equation 



dp 
dt 



+ V ■ {Ea) . (4.4) 



The term added to (3.5) to obtain in (4.4) is exactly the term which in [33] is 
added to the C-H equation to adapt it to a situation in which an external field is 
present. Problems and limitations to the use of macroscopic models like (4.4) are 
discussed in [19]: (4.4) in fact does not seem to be suitable to describe systems with 
strong (0(1) and not, like here, 0(7)) external fields. 

5. The hydro dynamic limit of the particle system: proofs 

The Proof of Theorem 1 is an immediate consequence of Theorem 2 and Propo- 
sition 1 given below. 

Some preliminary definitions:. In this section we will denote by e G Z*^ a 
unit lattice vector. By Ylie '^^ '^'^ mean the sum over all ej e Z'', {ej)i = 5j^i 
e {1, . . . , d}). Let us recall the Poisson rates 

c^{x, x + e;r]) = <^{(5/S.x,x+eH-y{r])) a: e A^, rj eil^ (5.1) 

where we used the notation 

for X, y e A^. We recall that $ e satisfies (2.9), $(0) = 1, and that c^(a;, y; ry) = 
if ja; — y| 7^ 1. In the proofs we shall denote by c^(x, y; rf) the rates in the case in 
which J = 0, that is 

c°(a;,a; + e;?7) = 1 a; e A^ (5.2) 

This process is called the simple exclusion process: the process with J ^ will be 
treated as a perturbation of the simple exclusion process. Moreover we will denote 
by c'P{x,y;r]) the rates when $(£') = exp(— i?/2) {E G R). The law of the simple 
exclusion process with initial condition fi is denoted by P^J'^ (E^''^), while the law 
of the process with rates c^(x, y; rj) is denoted by P^''^ (E^'^). 

Let Ml be the set of measurable functions p : T'^ ^ [0;1]- -^1 is equipped 
with the weak topology induced by duality by C{T'^), the continuous real- valued 
functions on T^, according to 



{p,G)= [ G(r)p(r)dr 
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where p e Mi and G e C{T^). Given 7 > we define the empirical measure of the 
process at r as 



where xa denotes the characteristic function of the set A C T^. The empirical 
averages of a function / : fi^ — > R over a ball of radius i? > is defined as 

^"'^'''f = \B(R) nZ<^\ ^ fi^^+yV) (5.4) 

where B{R) — {r : \r\ < R}, x E Aj and : Q-y — > Q-y is defined by {Tzr]){x) = 
'r]{x + z). For a e [0,1], let be the Bernoulli measure over such that 
v^{r]{x)) = Q! for all x e A^. We set 

/» = ^f(/)- (5.5) 

Given 77 G -D([0, 00), 0-y), the space of right continuous functions from [0, 00) to 
ri-y with limit from the left, equipped with the Skorohod topology, for every t G M"*" 
and r G T'^ we define 

p^(r,t) = p^(r;r]^-2t) (5.6) 

We will adopt the notation p^ — ■) and when we want to keep the time fixed we 
will write p^^t, which stands for p-y(-; 77-^-2^). So p^^t G Mi and py G D{[0, 00); Mi). 

With this notation, the family of measures {lJ''y}'y>o is an initial condition asso- 
ciated with pq if for any 5 > and any G C{T^) we have that 

hm E^;- (|(p-y,o,0) - (Po,0)| > 5) = 0. (5.7) 

(compare with (2.14)). 

A function p G C([0, cxd), Mi) is a weak solution of (2.16) if for all t G M+ and 
aU (f) G C2'i(T'^ X [0, 00)), p satisfies 



p{r,t)(j){r,t)dr — / po ('')</'('', 0)dr 
yd jj'd 



/ / p(r, s)ds(j){r, s)drds + / / p(?', s) A0(r, s)drds+ (5.8) 
J JQt J JQt 

1^ p(r, s)(l — p(r, s))(V J * p)(r, s)V(/)(r, s)drds 

J JQt 

in which Qt = T"^ x [0,t). 

Finally, for / G (/c G Z+), we set = ||/||co + E.„...,., Il^ii • ■■diJWco 

(zi, . . . , G and |zi + . . . + Zd| = A;), where || • \\co is the sup-norm. || • \\p will 
denote the L^-norm. 

11 



Theorem 2. For any t > 

lim P^;- (|(p^,„ 0) - (p„ 0)1 >5) = 0. (5.9) 

where pt is the unique solution of {5.8). 

Remark: a byproduct of the proof of Theorem 2 is that the random variable 
Py e -D([0, oo); Ml) converges weakly, as 7 approaches zero, to the deterministic 
trajectory p e C([0, 00); Mi), unique solution of (5.8). 

Proof of Theorem 2: For all e > 0, / : fi^ ^ R and r] e Q-y set 

ReAf'^v) = X]|Av,^-i,^(/;r7)-/(Av,^-i,^(7ro;r7)))| (5.10) 

X 

in which ttq : ^ {0, 1} is defined as 7ro(?7) = '?7(0). 

Let us recall the following result which is a straightforward extension to any 
dimension of Proposition 2.1 of [17]. 

Lemma 1. Given a cylindrical^ function f , t & M"*" and 5 > 



limsuplimsup7'*logP°''''''" ( 7'' / ReJf;r]s)ds > 5 \ = -00 (5.11) 

We have to extend Lemma 1 to the case of J 7^ and general initial condition. 
We are not going to extend the Lemma to the general case, but only to P^J''^^: 
in the general case we will obtain a weaker, but sufficient, estimate via a relative 
entropy argument. In order to extend Lemma 1, it is sufficient to show that there 
is a constant c such that 



I -^^i{Vt}te[o,,-^r]) I < c (5.12) 
for all r] e D([0, 7~^t]; fi-^), since (5.11) and (5.12) easily imply 

limsuplimsup7'^logPP'^^ (7'^ C R^Jf,r]s) > S] = -00 (5.13) 

e^O 7^0 \ Jo ' J 

under the same conditions stated in Lemma 1. The bound (5.13) will be obtained 
via direct evaluation of the Radon-Nikodym derivative, whose explicit expression 
is the following 

5^({^0tG[0,7-2r]) = 



■"^A function / on {0, +1}^ is said to be cylindrical if there exists a finite set A C 1,'^ such 
that fir]) = /(??') whenever ri(x) = rj' (x) for all x ^ A. Hence any cylindrical function has an 
obvious restriction to Jl-y, provided that 7 is small enough (i.e. if A C A^). The / which appears 
in formula (5.11) is precisely this restriction. 
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log I ^(^o) 1 - (5.14) 



/ y2Hix,x + e;r]t) -d{x,x + e;r]t)]dt+ (5.15) 



in which is the Poisson process that counts the exchanges of occupation 

number between x to x + e in the time interval [0, t]. We will bound the three terms 
of the RN derivative ((5.14), (5.15) and (5.16) ) separately. Let us observe that for 
(5.14), since {djjy/di'^^^){'r]) < 2'^ , the desired estimate is immediately obtained. 
Now and later on we will make use of the following Lemma. 

Lemma 2. For x E T'^ and rj E Q^, let 



K'^iv) = l^x,x+eH^{ri)-l{ri{x + e)-r]{x)) 



^''Y,v{z){e-VJ){-f{x-z)) 



(5.17) 



There exists a constant c\ such that for every a; e A^, every rj & and every e 

\h^^'%v)\ < ci^^ (5.18) 
in particular there is a constant ci such that for all x & Aj and all r] & 

\A^,x+eH^iv)\ < C27- (5.19) 

Remark: an immediate consequence of Lemma 2 is that the rates are bounded: 
for any J, /3 and $, there exists 70 such that 

sup sup Cy{x,y;rj) < 2. (5.20) 

7^(0,70) V'^^y 

Proof of Lemma 2: By definition 

-2A,,,i/^(r7) = Y.^rf'y{z)r^^'y{z') - r^{zUz')]^''j{{z - z')^) (5.21) 

z,z' 

add and subtract in the square brackets r]^'y{z)r]{z') and use the symmetry J(r) = 
J{—r) to get 

-2/^x,yH^{r^) = iviy) - r^{x)) Y.{v^'^{z) + r7(^))7''[J(7(^ - x)) - J{^{z - y))] 

z 

= -2{r]{y) - ?7(x)) J] vi^hViliv - ^)) - Jili^ - ^))] 

z (5.22) 

+ 2(r7(y) - v{x))V[J{0) - J{l{x - y))]. (5.23) 
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and if |a; — y| = 1 the modulus of the term (5.23) can be bounded by 2|| J||c7i7^"''^ 
and we can substitute the discrete gradient in (5.22) with (x — y) ■ V and the error 
will be bounded by || J||c'27^, so that (5.18) is proven with Ci = 2|| Jjlc*^ + || J||c2. 
Formula (5.19) follows immediately from (5.17) and (5.18). This ends the proof of 
Lemma 2. □ 

Let us go on with the proof of Theorem 2. The term in (5.15) can be written as 



I >J^*(^)(1 ~ + e)) [c^(x,a; + e;r7t) - l] dt 
By Lemma 2 we obtain that (5.24) is equal to 

\^Vnt{x){l - r]t{x + e)) 



(5.24) 



Jo 



x,e 



X 



+ r]t{x + e){l-r]t{x)) 



exp^-§7'+^5]r;,(z)(e-VJ)(7(^-a;-e)) + ^/.^+«'-(r7,)Vl } 



>dt 



(5.25) 



Set g^{x;r]) = {P/2){^'^^^'n{z){e ■ VJ){^{z - x)). The modulus of (5.25) can be 
bounded by 

- I - rit{x)rit{x + e) [exp(75r^(a;; rjt)) + exp(-75r^(a; + e; rjt)) - 2] 



+ r]t{x) [exp(75f^(a;;?7t)) + exp(-75f^(a;; ?7t))] \dt 
+ 2(ci72)(T7-2)(d7-^) <C7-'^ 



in which we used (5.20), (5.18), (5.19), the fact that {g-yix^r]) — g^{x + e;r/)| < 
(/?/2)|| J||c727 and c= c{d,T,P, \\J\\c^) < oo. The third and last piece (5.16) in the 
RN derivative is treated as follows. Given a sample path {r]t : t G [0,7~^t]}, call 
n = ^J^f^^^ the total number of exchanges in the interval [0,7~^t], which is 

finite for every r\. G -D([0, 7~^r], Mi). The i*^ jump (i = 1, . . . , n) happens at t = 
and it moves a particle from x^ to + e*. The last term in the RN derivative is 
equal to 



E 



a;,e 



log { + .,.+e ^ y^^ ,.+,.H^){v,-) 

^\c^^{x,x + e;rjt-)J ' ^ ^'^''i' 

=H^(r]t^-2) - Hj{r]o) 



(5.26) 



which is easily bounded by 2max^ \H^{r])\ < {\\J\\c°)l~'^ and so the proof of (5.12) 
(and so of (5.13)) is complete. 
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Clearly (5.13) implies that for e sufficiently small 

lim P^'/^- £ ' R,,^{f, r,s) >s\=0 



(5.27) 



The aim is to obtain a similar inequality for P^^ . To this purpose let us introduce 
the relative entropy between P^^^ and P^''^'' on the time interval [0,7~^t]: 



iy;(pf;^|pp'^^) = E^;^ 



log 



dp^;^ 
dp^ 



i{Vt}t&[0,-y-^r]) 



(5.28) 



and we will evaluate (5.28) directly by writing the RN derivative explicitely. The 
obtained expression is absolutely analogous to (5.15) and (5.16) (now the expres- 
sions — cP and c^/c^ are replacing — and c^/c^ respectively) and there is no 
term (5.14), since the initial condition is the same. Observe that by differentiating 
(2.9) we obtain 



(5.29) 



and so $'(0) = — ^»(0)/2 = —1/2. This implies that there exists a constant c such 
that \c^{x,y;r]) - cP{x,y;r])\ < c-f'^ and \\og{c^{x,y;r])/cP{x,y;r]))\ <c-f'^. Hence 



(P^;^|P^''^^) < ((i7"'^)(7"'T)(c72) + (c72)E^;^(n) < c'7 



(5.30) 



in which n is the total number of Poisson exchanges in [0,7~^t]. The bound in 
(5.30) follows directly from (5.20). By applying a well known entropy inequality 
(see e.g. [31], formula (2.18)) to our set up we obtain 



pM7 I .yd 



log 2 + 7^7 (P^;^|P^'''^) 



log 1 + l/P^^ (7^/0'" 'Re,M^Vs) > 5 



From (5.13), (5.30) and (5.31) directly follows that 



lim lim sup P^;-' r Re,M^Vs) >s]^0 

e^O \ Jo / 



(5.31) 



(5.32) 



€^0 

Let (p e C^'^(T'^ X [0, 00)): we will use the notation 4>t{-) = 4>{-, t). Define 

{P^,t, ^t)^ = 1'^^ V'r-Hix)(l){ix, t) (5.33) 



and observe that for any t > there exists c = c(</)) such that 



< 



Cf 



(5.34) 



which makes the two quantities appearing in the left-hand side of (5.34) inter- 
changeable for our purposes. Let us define 



15 



(5.35) 



where with partial derivative in t we mean the derivative with respect to the time 
dependence of (f), and 



9 2 

It is weU known (see for instance [7]) that 



(5.36) 



M, 



,(t; 0) = {pj^r, (f>T)j - (P7,o, 00)7 - / Ti^^{t)dt (5.37) 

^0 



and 



N^{T;cf>) = {M^{T;cf>)y 



(5.38) 



are P^'' —martingales with respect to the filtration generated by the family of ran- 
dom variables {?7t}t€K+- Moreover they vanish at r = 0. 

Let us start by computing Ti^y and r2,7. By using the definition of the generator, 
we obtain 

ri,7W = 

x,e 

+ ^''Y,M{lx,t)i^^-2t{x) (5.39) 



and the first term on the right-hand side of (5.39) can be rewritten as 

l'^Y.^A(j>){-lx,t)r)^-H{x)+ 



(5.40) 



X] ^7-2* - V-i-H^x + e)) ^ ^ r]^-2t{z){e ■ VJ){-f{x - z)) 



apart for corrections that are uniformly vanishing as 7 goes to zero and so jFi ,y(t)| is 
uniformly bounded (this follows from the observation following (5.29) that $'(0) = 
-1/2). Moreover 



r2,7W 

l~^^Cy{x,x + e;r]^-2t) 



x,e 
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E c^{x, x + e; r]^-2t) [(pi'yix + e), t) - (pi'yx, t)Y 



(5.41) 



x,e 



and from (5.41) we obtain 



|r2,7(t)|<2||</.||ci7'. 



(5.42) 



The fact that \{pj^ti 4't)\, |ri,7(t)| and |r2,7(t)| are bounded uniformly in t e [0,t] 
(r e M"*"), rj e D{[0,r];fl^) and 7 > 0, immediately implies that the family of 
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random variables {P'y}'y>o is (sequentially) relatively compact in L>([0, oo); Mi) 
(see for instance [17] and [15]). Moreover every limit point p of {P'y}'y>o lies in 
C([0, oo); Ml), which follows directly from the fact that with P^j^^ —probability one 
the sup^ I J^Jt'^^^ — J^-^^^) I = 1 and so the quantity in (5.33) can have jumps 
of size at most ||</>||co75 which vanish in the limit (see [7], Theorem 2.7.8). We are 
then left with the task of identifying the limit. 

Let us observe that (5.42) and Doob's martingale inequality imply that for every 
(5i > 0, e C^'i and r e M+ 



lim P^;-' sup |M^(t; 0)| > 5i = 0. 



7^0 



,te[o,r] 



(5.43) 



Observe that, because of the smoothness of and J, for every ^2 > 0, there is an 
eo > such that for all e e (0, eo] the following holds 



13 I ^'^^7ra;{vt)7^x+e{Vt){pj,t*e-VJ){^x){e-V(f)){^x,t)dt 



x,e 



~ I 7''y"Av,-i^,^(7ro7re,r/t)/5(p^,t *e- VJ)(7a;)(e- V0)dt 
and by (5.32) we obtain that for every 5 > 



<52 



lim lim sup V't'^ 



e-»0 



7-»0 



^ / 7'^y'Av,-i^,a;(7ro7re,r7t)(p^,t *e- VJ)(7a;)(e- V(^)(7a;,t)dt 



- / 7''^(Av,-i^,,(7ro,r7t))'[/3(p^,t*e-VJ)](e-V(/))(7x,Odt 



> 5 (5.44) 



<limlimsupP^;- /?||V0||co||J||ci 



X 



/ X] Ave-i^,a;(7ro7re, r]t) - (Av,-i^,^(7ro, r]t)f > 5 J =0 

•^0 a;,e ' ^ 



Take the limit in 7 along a convergent subsequence to obtain from (5.37), (5.40), 
(5.43) and (5.44) that for every 5 > with P^^'' —probability going to 1 as e goes to 
zero, the limit point p verifies the inequality 



yd 



p(r, T)0(r, T)dr — / p(r, O)0o('')dr — / / p(r, t)(9t0 + A0)drdt 

W4>(r, t){p - (tte * pf) [f3{p * V J)] drdt 



I Ic 



< 6 



in which = {2e)~'^X[-e,+e]''- e go to zero and by the arbitrariness of S we 
obtain that p is a weak solution of (2.16). To complete the proof of Proposition 1 
we have to show that there is only one weak solution. 
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Let us then consider pi and p2 weak solutions of (2.16) with the same initial 
datum. Set w = pi — p2- For each (ro,to) G Qt ^ind each e > let us define 

(l>ro,to,e{r^ t) = ^{to + e-t,r-ro) (5.45) 

in which 6'(r, t) — J^xeZ'^ Gtir — x) is the periodic heat kernel, which solves dt9 = 
on X (0,oo). So the test function (f)ro,to,e e C°°iQt,) solves 

dt(Pro,toAr, t) + A(/)„„t,„,(r, t) = (5.46) 

for all (r, t) G Qto- By using (5.46) we obtain (a"(p) = /3p(l — p)) 



/ 



I I ^Ao,toA^^ t) [(V J * w){r, t)a{p,{r, t)) + (V J * p2)(r, t)a'(p)«;(r, t)]drdt 
J JQto 

in which p e [(pi Ap2) VO, (pi Vp2) M]. Since a{pi) e [0,/3/4] and |c7'(p)| < /3 we 
obtain that 

/ (l)ro,toA'r^'to)w{r,to)dr 

<c[P\\VJ\\i] ess snp \w{r,t)\ [ [ |V(/),o,t„,e(r, t)|drdt (5.47) 

The integral term on the right hand side of (5.47) is bounded by c{^/to + e — ^/e), 
where c is a constant depending only on the dimension d. By observing that 
(pro,ta,e{^^ ^o) IS an approximate identity in e, we obtain that there exists a constant 
C (which depends only on the norm of V J and on d) such that for almost every 
(ro,to) £ Qt we have that 

\'w{ro,to)\ < C\/to ess sup \w{r,t)\ 

(r,t)eQto 

and so 

ess sup |t(;(r, t)| < C-\/t ess sup \w{r,t)\. 

{r,t)eQ-r (r,t)eQr 

Choosing r such that < 1 implies local uniqueness in time. But C depends 

only on J and d, so the global uniqueness follows by a bootstrap argument. □ 

Proposition 1. If po & C'^{T'^), then the solution p of (5.8) can he chosen in 
C2'i(T*^ X 



Proof: given a solution of (5.8) define 

F(r,t) = (VJ*pt)(r) (5.48) 

so F G (j<=c.,0(rpd ^ ■[]^+. j^d-j Observe that this means that p solves in the weak sense 
analogous to (5.8) the equation 

atp = V-{Vp-/3p(l-p)F} (5.49) 

with initial condition po G ■ But (5.49) is a nondegenerate parabolic equation 
which has a classical solution (Chapter 7, Section 4 of [11]), besides having a unique 
weak solution by the same argument used to prove uniqueness in the proof of 
Theorem 2. □ 
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